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Abstract 

It is shown that the modern equations of state for neutron star mat- 
ter based on microscopic calculations of symmetric and asymmetric 
nuclear matter are compatible with the lower bound on the maximum 
neutron-star mass for a certain range of hyperon coupling constants, 
which are constrained by the binding energies of hyperons in sym- 
metric nuclear matter. The hyperons are included by means of the 
relativistic Hartree- or Hartree-Fock approximation. The obtained 
couplings are also in satisfactory agreement with hypernuclei data in 
the relativistic Hartree scheme. Within the relativistic Hartree-Fock 
approximation hypernuclei have not been investigated so far. 



* Dedicated to Prof. Georg Sufimann on the occasion of his 70 birthday 
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Neutron star matter (NSM), bound by gravity, differs from high density nu- 
clear matter produced in heavy ion reactions in several respects: 1) Since the 
repulsive Coulomb force is much stronger than the gravitational attraction, 
NSM is much more asymmetric than terrestrial matter in heavy ion colli- 
sions. 2) The weak interaction time scale is small in comparison with the 
lifetime of the neutron star (NS), but large in comparison with the lifetime 
of high-density matter in heavy ion reactions. For that reasons dense matter 
in high-energy reactions has to obey the constraints of isospin symmetry and 
strangeness conservation whereas NSM is a charge neutral system with no 
strangeness conservation (see, for instance, Ref. [||] with further references). 
Extraterrestrial NSM is for that reasons a rather theoretical object with a 
very complex structure. Since its density stretches over an enormous range, 
reaching from crystalline iron on the surface to several times of nuclear matter 
saturation density in the core, one encounters in the theoretical descriptions 
many obstacles, for which one has to introduce theoretical assumptions and 
extrapolations [1], [|. However one can try to combine both high density 
systems in a common theory, which is possible in modern field theoretical 
approaches. The first systematic investigation in this respect was performed 
by Glendenning who used the standard nuclear relativistic Lagrangian 
extended by inclusion of hyperons and deltas, i.e. 

C(x) = E C °b{x) 

B=p,n,E±°,A,S° -.,A->°.+>++ 

r i (1) 
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where the interaction is mediated by the mesons a, u, ■ ■ • and the leptons 
are treated as free particles. Furthermore a scalar meson self-interaction is 
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included. More explicitly the hadron part is given by: 
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with 



F Z»(x) = d^ix) - dvWpix) , F p ^(x) = d^p v (x) - d v p v {x) . (3) 

Within this scheme the vast majority of investigations was performed using 
coupling parameters in the nucleonic sector adjusted to the nuclear matter 
parameter (p , E/A, K, J) and the Dirac mass rh and a variety of choices 
for the hyperon couplings (see, for instance, Refs. |], |3|, [|, ||). The calcu- 
lations for such Lagrangians could be performed without numerical conver- 
gence problems and resulted in a high strangeness contents of NSM |], |3|. 
However it was pointed out in two more recent investigations, that due 
to the selected high Dirac mass (~ 0.79 m^) such parametrizations are 
neither in accordance with standard phenomenological parametrizations ||, 
as, for instance, the frequently employed parametrizations NL1 and NL-SH, 
or with parametrizations adjusted to the outcome of microscopic relativis- 
tic Brueckner-Hartree-Fock calculations (RBHF), respectively, which give a 
lower Dirac mass of approximately 0.6 0? @] (necessary also, for instance, 
for a correct reproduction of the spin-orbit-splitting 0). Attempts to use 
directly such parametrizations in NSM-calculations, which one should prefer 
theoretically, have not been successful, since they lead to negative nucleon 
Dirac masses by inclusion of hyperons 0, 0. This strong decrease of rh, 
caused by the stronger ex-field in the case of lower Dirac masses at satura- 
tion in comparison to the high choice of rh by Glendenning, seems to be also 
not in accordance with the behaviour of rh in RBHF-calculations for higher 
densities in symmetric nuclear matter ||. For that reason one has to correct 
this feature of the RH, which occurs also in the relativistic Hartree-Fock- 
approximation, by a slight change of the Lagrangian for higher densities, 
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which gives a more moderate decrease in this region similiar to the RBHF 
values but does not change the equation of state (EOS) in the region up to 
moderate densities above saturation nuclear density PJ. Within this scheme, 
described in detail in Ref. JD], where the changes in the masses are moderate 
and the high-density behaviour is through the dominance of the cu-repulsion 
hardly altered, one obtains now in the RH-approximation NS-matter com- 
positions with a high strangeness component which are qualitatively simi- 
lar to former investigations, but have now the advantage, that they are 
based either on a correct phenomenological description of nuclear proper- 
ties and the outcome of microscopic RBHF-calculations, with an additional 
assumption, done in a controlled and changeable manner by one additional 
parameter. A closer inspection reveals that all the other previous attempts, 
where such problems did not arise at a first glance, have made more or less 
implicitly, for instance, by selecting a high Diracmass || [5| or by introduc- 
ing a io self- interaction ||10|| , such an additional assumption, however with a 



rather specific enforcement of the high density behaviour on the Dirac mass 
(for details, see Ref. [![]). In the RH-scheme, A's usually do not occur, since 
due to the large p-coupling, needed to reproduce the correct symmetry coef- 
ficients, the charge-favoured but isospin-unfavoured A~ is suppressed [[I], |3|. 
Since the data from the hypernuclei are not yet conclusive with respect to the 
hyperon couplings (for more details, see Refs. [D], [12]]), it is common in more 



modern investigations to use the SU(6) symmetry for the vector couplings 
and adjust the a-coupling with respect to the lowest hyperon level in nuclear 
matter. However the data from hypernuclei still permit a large bandwidth 
for the couplings |], H 01 [H> HUl- As it was pointed out by Glendenning 



and Moszkowski |13| one can also use the maximum neutron-star mass of 
at least 1.5 solar mass, the so-called Oppenheimer-Volkoff limit Mov, as 
an additional constraint besides the binding energy of the hyperon in mat- 
ter to fix the permitted pairs of coupling constants x a = g a H / 9un{< 1) and 
%u = 9uiH I '5Wv(— !)• However one has to keep in mind that such fixations de- 
pend strongly on the chosen nucleonic parametrization in combination with 
the selected approximation PJ. It seems therefore worthwhile to revisit this 
problem in a short note in view of the discussed theoretically more satisfy- 
ing parametrizations, based on the outcome of RHBH-calculations with the 
modern Brockmann-Machleidt OBE-potentials QTJ, |7|, [14]]. The parametriza- 
tions are given in Table I and Table II. 



More subtle is the case of the RHF-approximation H, 14], where the con 



straint of the hyperon binding energy in nuclear matter is difficult to imple- 



4 



ment [I]]. A further new feature is that due to the exchange contributions 
the g p -coupling becomes smaller in the direction of the values of the realistic 
OBE-potentials with the net effect that now the A's can play an impor- 
tant role in the composition of NSM, so decreasing the strangeness contents 
0. The A-component is rather sensitive to the choice of the A-coupling, 
and one obtains qualitatively rather different compositions, reaching from 
the standard picture for lower A-couplings, as recommended by Rapp et al. 
[ 151 ) to A-dominated compositions with much smaller strangeness obtained 
f° r 9a/ 9n = 10- Since a final solid statement about the choice of these 
couplings cannot be made, one cannot exclude, at present, the possibility 
that the composition of NSM might deviate from the standard picture of 
a high strangeness component, which may be a special feature of the RH 
approximation with a large p-coupling (for details, see, Ref. 0). The chosen 
parametrizations are given in Table III. The dependence of the relative scalar 
coupling for the hyperons for a selected relative vector coupling is exhibited 
in Fig. 1. As far as the RH-approximation is concerned all parametrizations 
with x w > 0.5 give maximum masses above 1.6 M , the lowest maximum 
gravitational masses are obtained for no coupling, i.e. RH6, where one ob- 
tains Mqv — 1-17 Mqv- The Oppenheimer-Volkoff masses for the SU(6)- 
values are: Mov — 1-68 M & with a baryonic mass of Movb = 1-92 M Q 
for nonrotating NSs and M ov = 2.03 M g (M OV b = 2.32 M ) for NSs ro- 
tating with their Kepler frequency Qk — 7.97 x 10 3 s -1 . The OV-masses 
increase with stronger hyperon-coupling to approximately 2.15 M Q (Q = 0) 
or 2.45 M & (Q = Qk), respectively, for the EOS-RH4. In this case with the 
strongest hyperon interaction one obtains the largest hyperon contribution, 
which is shown in Fig. 2. Since the masses of most of the pulsars are in 
the range of 1.5 M & (the most accurately measured mass is that of the PSR 
1913+16 with M/M Q = 1.422±0.003 |||, all pairs of hyperon couplings seem 
to be compatible or at least not in contradiction with the observations. The 
comparison for the RH is depicted in Fig. 3. In the framework of the RHF- 
approximation |], [14] the situation is more complicated. Naive use of the 
hyperon couplings from the RH-treatment is not compatible with the binding 
energies of hyperons in saturated nuclear matter, since in this case the ex- 
change contributions contribute strongly different in the time-like and scalar 
parts of the self-energies and hence cause an overestimation of the A-binding 
energy in matter. In this case the choice with SU(6) vector couplings is not 
sufficient to reach a OV-mass of 1.5 M & (Mqv — 1-35 M & ). By correcting 
this feature one obtains lower scalar couplings (see Table III) in comparison 
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with the RH. The OV-masses corresponding to the SU(6) A-coupling, i.e. 
RHF3, give now M ov = 1-59 M & {M OVB = 1.86 M ) for a static, nonro- 
tating spherical NS and M ov = 1-88 (M OV b = 2.20 Af ) for the deformed 
NS rotating at Kepler frequency Qk — 10.51 x 10 3 s _1 . For this EOS the As 
contribute significantly, since the strong coupling of the As in combination 
with the weaker p-coupling is capable to overcome their iso-spin unfavoured 
deficit. However decrease of the A-coupling according to Ref. fl5j is sufficient 



to push out again the iso-spin unfavoured As and the hyperons contribute in 
a similiar manner as in the RH. The results for the OV-masses (RHF9) are 
then in the same range, namely: Mov — 1-65 M (Movb = 1-89 M ) for the 
static NS and M ov = 1-94 M & (M OV b = 2.22 M ) for the rotating star with 
the Kepler frequency Qk = 8.26 x 10 3 s _1 . The outcome for the RHF is de- 
picted in Fig. 3, too. From figures 1 and 3 one can infer that the parametriza- 
tions with 0.5 < x w \ < 0.8 for the RH- and RHF-approximations are with 
respect to the hyperon couplings in accordance with the NS-mass data, since 
they are capable to sustain a mass above 1.5 M , which is in accordance with 
the observation of the pulsar PSR 1913-16 with a mass of 1.44 M . In this 
context one should remark that for an isolated NS the maximum mass is 
smaller than the OV-mass emerging from the calculation of a cold NS. This 
difference is caused by the fact that NSs at birth are composed of so-called 
supernova matter with a high lepton fraction characterized by adjusted to 
the outcome of our RBHF-calculations one needs a slightly lower x a \ = 0.59 
in order to obtain the correct lowest binding energy of A in 17 O, PD| . In 
view of the included density-dependence such deviations seem reasonable 
| 20| . For the RHF-approximation a calculation of A-hypernuclei has not 



been performed so far. 

In conclusion we have determined a range of possible EOSs for NSM, based 
on microscopic EOSs of symmetric and asymmetric nuclear matter, not pos- 
sessing the unjustified high Dirac masses of earlier treatments and which are 
mutually compatible with neutron-star masses, A-binding in nuclear matter 
and hypernuclear levels (in the RH). The relation between the scalar and 
vector coupling was fixed by the binding energy, for instance, of the lowest 
A-level in nuclear matter by the Hugenholtz-Van Hove theorem. So far the 
comparison has been only tested for the relativistic Hartree-approximation, 
since hypernuclei calculations in the RHF-approximation are not available 
at present. 
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TABLE I 





9a 


9u 


9 P 


10 3 x b N 


10 3 x c N 


fp/g P 


RHA 


9.58096 


10.67698 


3.81003 


3.333665 


-3.52365 




RHFA1 


9.28353 


8.37378 


2.10082 


3.333689 


-2.15239 




RHF A2 


9.24268 


8.25548 


2.19809 


2.96514 


-2.68614 


3.7 


RHF A3 


9.16665 


8.07540 


2.37987 


1.95524 


-2.36335 


6.6 


RHB 


9.59169 


10.68084 


3.66541 


3.62616 


-4.17140 




RHFB1 


9.36839 


8.40466 


1.77326 


3.74354 


-3.18456 




RHFB2 


9.33266 


8.32154 


1.86078 


3.44306 


-3.46261 


3.7 


RHFB3 


9.26782 


8.19391 


2.02216 


2.67292 


-3.18198 


6.6 



TABLE II 



Approximation 










X<tA 


x wA 


RH1(SU(6)) 


0.611 


2/3 


0.346 


1/3 


1 


1 


RH2 


0.40 


0.40 


0.40 


0.40 


1 


1 


RH3 


0.64 


0.7 


0.63 


0.7 


1 


1 


RH4 


0.79 


0.9 


0.79 


0.9 


1 


1 


RH5 


1 


1 


1 


1 


1 


1 


RH6 





















TABLE III 



Approximation 


XaZA 








%aA 


XujA 


RHF1 


0.611 


2/3 


0.346 


1/3 


1 


1 


RHF2 


0.79 


0.9 


0.79 


0.9 


1 


1 


RHF3 


0.44 


2/3 


0.26 


1/3 


1 


1 


RHF4 


0.36 


0.5 


0.35 


0.50 


1 


1 


RHF5 


0.46 


0.7 


0.45 


0.7 


1 


1 


RHF6 


0.56 


0.9 


0.55 


0.9 


1 


1 


RHF7 


1 


1 


1 


1 


1 


1 


RHF8 


0.44 


2/3 


0.26 


1/3 


0.625 


0.625 


RHF9 


0.44 


2/3 


0.26 


1/3 


0.625 


1 


RHF10 


0.44 


2/3 


0.26 


1/3 


0.75 


1 


RHF 11 


0.44 


2/3 


0.26 


1/3 
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Table captions 



Table I: Parametrizations of the RH and RHF-Lagrangian adjusted to the 
RBHF-calculations. For the masses the following values were selected 
(MeV): rriN = 939, ra a = 550, m w = 738, m n = 138, m p = 770 (g n = 
1.00265, f%/4:7t = 0.08). The parametrizations are labelled as follows: 
RHA = RBHA etc. for the potential A (effective mass at the Fermi 
surface m = 617.8 MeV (A); 621.8 MeV (£)). 

Table II: Relative coupling strengths of the hyperons in the relativistic Hartree 
approximation. In Ref. |1] occured a misprint for x a ^(= 0.676) and 
x a s(— 0.342). For RH1 - RH4 x a u was adjusted for a selected x^h 
according to the Van-Hove-Hugenholtz theorem, which gives for the 
hyperon potential depth U(H) = x a n^s + Xu>h^o, where S s and S 
are the scalar and vector part of the self-energy at saturation. For 
comparison we included also the universal coupling (RH5) and the free 
Fermi gas approximation for the hyperons (RH6). For the nucleonic 
sector the parametrization RHB was used. 

Table III: Relative coupling strengths of the hyperons in the relativistic Hartree- 
Fock approximation. In RHF1 and RHF2 the coupling strengths of the 
RH-approximation were used. For RHF3 - RHF11 x a was adjusted for 
a given x u according to the binding energy of the hyperon in symmetric 
nuclear matter at saturation. For the nucleonic sector the parametriza- 
tion RHF Bl was used. RHF10 is characterized by fixing the potential 
depth of A - to the same value as in RH1 (for the composition, see, 
Ref.!). 

Figure captions 

Fig.l. Constraint of the relative a-coupling constant for a given cu-coupling 
due to the lowest hyperon binding energy in saturated nuclear matter 
for the RH- and RHF-approximation, respectively. 

Fig. 2. Baryon-lepton composition as function of the density of a NS in the 
RH-approximation with the strongest hyperon-coupling (RH4). 

Fig. 3. Oppenheimer-Volkoff masses for static neutron stars for the different 
approximations as function of the relative a;^^— coupling. X(T h is ad- 
justed for these parametrization according to the lowest hyperon bind- 
ing energy in matter. The horizontal line gives the mass of the pulsar 
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PSR 1913+16. x^h = 0.8 (x aH < 0.72) gives the highest value for the 
coupling constant compatible with properties of hypernuclei [|T^, [ZD] ]. 
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